% this function is used for merging the 'r0m(m=1,2,...)' into 'r0'
function r0 = r0m_merge(r0m)
% r0m: r0 for each layer of turbulence path
% r0: the total Fried constant

r0 = (sum((r0m.^(-5/3))))^(-3/5);

% end of the function
end
